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Abstract — The greatest changes in elevation in Greenland and 
Antarctica are happening along the margins of the ice sheets where 
the surface frequently has significant slopes. For this reason, the 
upcoming Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) 
mission utilizes pairs of laser altimeter beams that are perpendic- 
ular to the flight direction in order to extract slope information in 
addition to elevation. The Multiple Altimeter Beam Experimental 
Lidar (MABEL) is a high-altitude airborne laser altimeter de- 
signed as a simulator for ICESat-2. The MABEL design uses mul- 
tiple beams at fixed angles and allows for local slope determina- 
tion. Here, we present local slopes as determined by MABEL and 
compare them to those determined by the Airborne Topographic 
Mapper (ATM) over the same flight lines in Greenland. We make 
these comparisons with consideration for the planned ICESat-2 
beam geometry. Results indicate that the mean slope residuals 
between MABEL and ATM remain small « 0.05°) through a 
wide range of localized slopes using ICESat-2 beam geometry. 
Furthermore, when MABEL data are subsampled by a factor of 
4 to mimic the planned ICESat-2 transmit-energy configuration, 
the results are indistinguishable from the full-data-rate analysis. 
Results from MABEL suggest that ICESat-2 beam geometry and 
transmit-energy configuration are appropriate for the determina- 
tion of slope on ~90-m spatial scales, a measurement that will be 
fundamental to deconvolving the effects of surface slope from the 
ice-sheet surface change derived from ICESat-2. 

Index Terms — Airborne lidar, laser altimetry, photon counting. 

I. Introduction 

T HE Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) 
is a NASA mission scheduled to launch in 2016. ICESat- 
2 is a follow-on mission to ICESat (2003-2009) and will 
extend the time series of elevation-change measurements aimed 
at estimating the contribution of polar ice sheets to sea-level 
rise, quantifying changes in sea ice freeboard, and enabling 
determination of global vegetation height. Relative to ICESat, 
ICESat-2 represents an opportunity to improve estimates of ice- 
sheet surface-elevation change by resolving surface slope. 

ICESat carried the Geoscience Laser Altimeter System 
(GLAS), a single-beam laser altimeter that recorded the re- 
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Fig. 1. (Left) Schematic ICESat, (center) ICESat-2, and (right) MABEL 
(dashed lines) beam geometry and (gray lines along ice-sheet surface) reference 
ground tracks. MABEL allows for beam-geometry changes with a maximum 
ground spacing of d= 1 km at 20 km. 

ceived laser energy as a waveform to determine surface eleva- 
tion, based on reflected 1064-nm wavelength laser light with 
a 40-Hz pulse-repetition frequency. GLAS sampled 50-70-m- 
diameter footprints every ^172 m along a series of repeated 
tracks [1], providing information about surface elevation and lo- 
cal slope along track but not across track. ICESat ground tracks 
did not repeat exactly, further complicating the GLAS determi- 
nation of surface-elevation change; the postprocessed footprints 
were offset from a reference ground track by up to 200 m. Thus, 
for a single-beam laser altimeter, such as ICESat, change in 
surface elevation with respect to time is difficult to separate 
from the effects of local surface slope. Ultimately, to make this 
separation, ICESat required up to 12 passes in areas with either 
low slope or subtle elevation change, such as East Antarctica. 

ICESat-2 will carry the Advanced Topographic Laser Al- 
timeter System (ATLAS), which uses a different detection 
strategy than GLAS. Specifically, ATLAS will be a six-beam 
photon-counting laser altimeter using short (< 2 ns) 532-nm 
wavelength pulses, with a 10-kHz repetition rate. In a photon- 
counting system, single-photon sensitive detectors are used, and 
the arrival time of any detected photon is recorded. ATLAS 
will have a ^10-m-diameter footprint and a ^70-cm along- 
track sampling interval [2]. An accurate assessment of ice-sheet 
elevation change is dependent upon knowledge of local slope 
[3]. Therefore, the six beams are arranged in three sets of pairs. 
Spacing between the three sets of pairs is ^3 km to increase 
sampling density, whereas spacing between each beam within 
a given pair will be ~90 m to make the critical determination 
of local slope on each pass (see Fig. 1). Therefore, elevation 
change can be determined from only two passes of a given area. 
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Fig. 2. Location map of the 50-km along-track Summit Area site and the 
150-km along-track Southern Traverse site on the MODIS Mosaic of Greenland 
[17]. (Inset) Operation IceBridge P-3B captured in the NASA ER-2 Cirrus 
Digital Camera System [18] at the black 

Given this new approach to spacebome surface elevation 
measurement, an airborne instrument, the Multiple Altime- 
ter Beam Experimental Lidar (MABEL), was developed to: 
1) enable the development of ICESat-2 geophysical algorithms 
prior to launch; 2) provide detailed error analysis of the ATLAS 
measurement strategy; and 3) provide ATLAS model valida- 
tion. MABEL is a photon-counting multibeam lidar sampling 
at both 532- and 1064-nm wavelengths using short (^1.5-ns) 
laser pulses. MABEL beams are arranged in a linear array, 
perpendicular to the direction of flight. The system allows for 
beam-geometry changes between flights with a maximum view 
angle of ±1 km from a 20-km nominal altitude achieved during 
the 2010-2012 deployments using a NASA ER-2 aircraft (see 
Lig. 1). The repetition rate of MABEL is variable (between 5 
and 25 kHz); most flights during the 2010-2012 deployments 
used 5 kHz. At this nominal altitude, repetition rate, and an 
aircraft speed of ^200 m/s, MABEL samples a ^2-m footprint 
every cm along track [4]. During these initial MABEL 
deployments, beam geometry (specifically the spacing between 
the individual beams) was configured to mimic ICESat-2 [5]. 
Following engineering test flights in December 2010 and March 
2011, MABEL was deployed to Greenland in April 2012 to 
collect data over polar targets (see Fig. 2). 

Operation IceBridge is a NASA airborne campaign intended 
to bridge the data gap between ICESat and ICESat-2. Opera- 
tion IceBridge hosts a suite of instruments [6], including the 
Airborne Topographic Mapper (ATM). ATM is a lidar that 


conically scans at a rate of 20 Hz, with an off-nadir scanning 
angle of ^15°. Like GLAS, ATM digitizes returned energy as 
a waveform with derived surface elevations based on 532-nm 
wavelength pulses and a 5-kHz pulse-repetition frequency [7], 
[8]. The ATM flights used in this letter were conducted using 
NASA P-3B at an aircraft speed of ^100 m/s, with a nominal 
elevation of 500 m above ground level. At this air speed, eleva- 
tion, and repetition frequency, ATM generates a 1-m footprint 
and a scanning swath width of ^250 m. 

Logistics and cloud-free weather allowed for coordinated 
surveys between ATM and MABEL over the Greenland Ice 
Sheet (see Fig. 2). Here, we present MABEL multibeam de- 
termination of the ice- sheet surface and compare it with that 
determined by ATM, including local slope assessments. These 
comparisons are made with consideration for the ICESat-2 
planned beam geometry and relative signal strength. 

II. Data and Methods 

Both MABEL and ATM simultaneously surveyed a 150-km 
“Southern Traverse” of the Greenland Ice Sheet on April 20, 
2012 (see Fig. 2). Additionally, MABEL made three passes 
over a 50-km stretch of ICESat track 0412 in the vicinity of 
Summit Station, Greenland, on April 8, 2012. ATM made a pass 
of the same ground segment on April 11, 2012 (see Fig. 2). 
This ground segment has been used as a calibration site for 
ICESat [9]. 

The NASA IceBridge ATM Level-2 Icessn Elevation, Slope, 
and Roughness (ILATM2) for April 11 and 20, 2012, were 
obtained from the National Snow and Ice Data Center (NSIDC; 
http://nsidc.org/data/ilatm2.html). This is a resampled and 
smoothed elevation data set [10] that provides four across- 
track elevations per timestamp every ^35-m along track, which 
allowed for the trivial calculation of across-track slope. The 
total across-track span for this data set, for the flights used in 
this analysis, was approximately 150 m. 

MABEL data (release 8) for April 8 and 20, 2012, were 
obtained from the NASA ICESat-2 website http://icesat.gsfc. 
nasa.gov/icesat2/). Each data file contains 1 min of data for 
every available beam. The data files contain photon arrival 
times resulting from reflected laser light (i.e., signal photons) 
and background photons due to sunlight (i.e., noise photons). 

In order to discriminate coarse signal photons from noise 
photons and derive ice- sheet surface elevation, we developed 
an algorithm based on histograms of the photon data. Evolving 
from techniques applied to other photon-counting lidars, such 
as Slope Imaging Multi-polarization Photon-counting Lidar 
[11], we generated 0.125-s (^25 m; 625 shots) along-track 
segments of data and histogrammed the photon data at 10-m 
vertical resolution. We identify signal photons in four sequen- 
tial steps. First, the mean and standard deviation of total bin 
counts are calculated. Bins with more counts than the mean 
plus the standard deviation are considered to potentially con- 
tain signal photons. Second, we exclude the potential signal 
bins found in step 1 and recalculate the mean and standard 
deviation for the remaining bins. These values primarily reflect 
the mean and standard deviation of the background count rate. 
We add bins with more than 15 times the standard deviation 
of background counts to the potential signal bins found in 
step 1 . Third, we include any individual photon that is within 
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±10 m of signal photons found in steps 1 and 2 as additional 
signal photons. Last, we select the median photon elevation, 
plus all photons within ±0.5 m of the median elevation, as 
signal photons. This approach assumes that the background 
photons are randomly distributed and that the returned pulse is 
symmetric. Data from surfaces with large slopes or roughness 
can generate a returned pulse wider than the 1-m window 
of step 4 (~3.3 ns of flight time, or about 2.5 times the 
transmitted pulsewidth). Additional analysis of MABEL data 
from the ice-sheet interior indicates that the returned pulse 
shape is largely symmetric and less than twice the transmitted 
pulsewidth [12]. 

For direct comparison of the surveys, ATM tracks and 
MABEL beams were chosen to most closely mimic the 90-m 
spacing of the ICESat-2 beam geometry. ATM tracks 2 and 3 
were used for this analysis because they have ground separation 
of ~85 m. Elevations from tracks 2 and 3 from the same along- 
track time were then used to calculate ATM across-track slope. 
MABEL beam 6 (center of the array) and beam 5 (^85-m 
ground spacing from the center of the array) were chosen for 
analysis as they have an across-track ground separation similar 
to the ATM tracks used. To determine MABEL across-track 
slope, signal photons from beams 5 and 6 were interpolated 
along track to a common time so that, similar to ATM, an 
across-track slope could be then calculated for each increment 
of along-track time. The across-track slopes for both ATM and 
MABEL were then compared. 

MABEL beams have variable signal strengths; however, 
beams 5 and 6 are the most similar to the expected radiome- 
try of the strong beams of ATLAS. Along-track data density 
differed within and between flights based on variables that 
affect reflectivity, including weather conditions, time of day, 
and sun-incidence angle. For the data used in this analysis, 
the full-rate along-track data density average for both beams 
was always greater than 4 signal photons/m. For the Southern 
Traverse flight, the along-track data densities were 3.4 and 3.9 
signal photons per 70 cm for beams 5 and 6, respectively. For 
the Summit Area flight, data densities were 3.1 and 3.4 signal 
photons per 70 cm for beams 5 and 6, respectively. 

A strong-beam/weak-beam pair will be used for ICESat-2 
slope determination; the energy associated with the weak beam 
will be reduced by a factor of 4. Therefore, the expected number 
of signal photons per laser shot (every 70 cm along track) 
between the strong beam and the weak beam will also differ 
by a factor of 4. The current best estimates of expected signal 
photons per laser shot vary with season and surface type. Based 
on ICESat-2 engineering models, under similar conditions as 
the 2012 MABEL survey, we expect ICESat-2 to record 8.5 
and 2.1 signal photons every shot (or 70 cm along track) for the 
strong and weak beams, respectively [13]. Thus, the MABEL 
full-rate data used in this analysis suggest data densities of 46% 
of the expected ICESat-2 data densities. MABEL engineers are 
currently working to increase signal strength to achieve the 
expected ICESat-2 data densities, which will facilitate more 
direct MABEL to ATLAS comparisons. 

To further assess accurate ground characterization given the 
ICESat-2 planned configuration, all photons associated with 
one of the MABEL beams (beam 5) were subsampled by a 
factor of 4 and then reprocessed through the ground- finding 
algorithm, to simulate the expected radiometric relationship be- 


tween the ATLAS strong and weak beams. After subsampling, 
data densities were 0.9 and 0.8 signal photons per 70 cm for the 
Southern Traverse and Summit Area, respectively. To determine 
MABEL across-track slope, the ground-signal photons from 
beam 6 and the subsampled ground-signal photons from beam 
5 were again interpolated to a common time so that an across- 
track slope could be calculated as described above. Therefore, 
the beam with the fewest along-track samples (the weak beam, 
5) limited the total number of samples that was used in the slope 
determination. 


III. Results 

A. MABEL Accuracy and Precision 

ATM is calibrated each flight by comparing airborne lidar 
elevations with elevations derived from a ground-based GPS 
survey of a portion of the departure airport apron (in this case, 
a 1.4-km along-track portion of the Kangerlussuaq, Greenland 
airport apron). The Kangerlussuaq GPS ground survey data 
(differentially postprocessed with an RMS < 5 cm) were pro- 
vided by the ATM project [14], which allowed us to make an 
assessment of the ATM precision for flights on April 1 1 and 20, 
2012 (± 16.3 cm). MABEL overflew the apron during the April 
20, 2012, flight, allowing for the determination of the accuracy 
and precision of the elevations derived from the ground-finding 
algorithm for both of the 532-nm beams used in this analysis 
(beams 5 and 6: —1.99 m ± 14 cm and —1.76 ± 14 cm, 
respectively). MABEL elevations from each beam were then 
corrected independently based on these biases. 

B. MABEL Local Slope Assessment 

Along the Southern Traverse, ATM and MABEL both col- 
lected data along a transect across the Greenland Ice Sheet, 
including a wide range of elevations and slopes (see Fig. 3); 
in general, the eastern flank of this traverse is much steeper 
than the western flank. The ground surface detected by MABEL 
traced that determined by ATM [see Fig. 3(b)]. Visual assess- 
ment of the across-track slope results indicated that relative to 
ATM, MABEL is able to determine the character of the surface 
throughout a wide range of slopes [see Fig. 3(c)]. The mean 
slope values along the Southern Traverse for ATM and MABEL 
were 0.99° and 1.03°, respectively, with a mean slope residual 
between the two sensors of 0.03°. 

Similar results were observed at the Summit Area site, which 
is very flat relative to the Southern Traverse flight line (see 
Fig. 4; note the change in the scale of the 7/-axes, relative to 
Fig. 3). Again, the ground surface detected by MABEL traces 
that determined by ATM [see Fig. 4(a)]. The mean slope values 
for the Summit Area for ATM and MABEL were 0.07° and 
0.1 1°, respectively, with a mean slope residual between the two 
sensors of 0.04°. 


C. Elevation Uncertainty 

Surface elevation uncertainty is related to geolocation (hori- 
zontal) uncertainty and local slope angle a [15] 

elevation uncertainty = tan(<r) x horizontal uncertainty. (1) 
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Fig. 3. ATM and MABEL data, from the same day, along the Southern 
Traverse sampling site (see Fig. 2). (a) MODIS Mosaic of Greenland [17] with 
(red) ATM and (blue) MABEL flight paths, (b) ATM and MABEL elevations 
along the flight paths, (c) ATM and MABEL across-track slope angles in 
degrees. 


Fig. 4. ATM and MABEL data along the Summit Area sampling site (see 
Fig. 2). (a) (red) ATM and (blue) MABEL elevations along the ground tracks, 
(b) ATM and MABEL across-track slope angles in degrees; note that the scales 
of the y-axes are substantially different from Fig. 3. 

MABEL has geolocation knowledge of 0.1 mrad (0.0057°) 
[16]; for a 20-km nominal altitude, this produces a 2-m hori- 
zontal uncertainty on the surface. Based on (1) and the range of 


Fig. 5. Across-track slope angles based on ATM and simulated weak-beam 
MABEL along the (a) Southern Traverse and (b) Summit Area sampling sites. 
Across-track slope angles for (red) ATM and MABEL subsampled by factors 
of (blue) 4 and (black) 10; at a subsampling rate of 10, the standard deviation 
of the slope differs. 

slope values for the Summit Area [with a minimum of ~0.1° 
from Fig. 4(b)] and the Southern Traverse [with a maximum 
of 6° from Fig. 3(c)], MABEL geolocation-induced elevation 
uncertainties for this letter will range from 0.4 to 21 cm, 
respectively. 

The other major sources of elevation uncertainty are due 
to errors in the surface-finding algorithm and time-tag errors 
within MABEL. Laboratory evaluation shows that the tim- 
ing precision of MABEL is ^200 ps (or 3 cm of range). 
While the contribution from the algorithm is difficult to defini- 
tively assess, it is significantly less than the return pulsewidth 
(^1.8-ns full- width at half maximum); therefore, we assume a 
value of 15 cm. Taking the square root of the sum of the squares 
of these terms leads to a cumulative elevation uncertainty range 
of 15-26 cm. 

D. MABEL Subsampled Assessment of Local Slope 

The slopes based on the MABEL ground signal subsampled 
by a factor of 4 were robust with respect to ATM slope determi- 
nation and the full-rate MABEL ground signal [compare Fig. 5 
to Figs. 3(c) and 4(b)]. The across-track slopes derived from 
subsampled ground-signal photons, for both the variable terrain 
of the Southern Traverse and the flat terrain of the Summit Area, 
were statistically indistinguishable from the MABEL full-rate 
across-track slopes. 

When the photons associated with full-rate beam 5 were 
subsampled by a factor of 10, substantially reducing the data 
rates of the simulated weak beam (to 0.3 signal photons per 
70 cm for both the Southern Traverse and the Summit Area), 
relative to ATM and the full-rate MABEL data, the mean slope 
values for both study areas remained relatively unchanged (see 
Fig. 5). However, the standard deviations of the slopes differed. 
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For the Southern Traverse, the mean slope slightly decreased 
from 1.03° to 0.98°, whereas the standard deviation changed 
from 1.14° to 2.41°. For the Summit Area, the results were 
more dramatic; the mean slope decreased from 0.1 1° to —0.43°, 
whereas the standard deviation changed from 0.09° to 5.87°. 
These changes in the standard deviation of the data (seen as 
outliers in Fig. 5) suggest that the subsampling of the full- 
rate MABEL photon data is having an adverse effect on the 
ground- finding algorithm. While the ground- finding algorithm 
could be modified to search for the surface on different along- 
track spatial scales, these results give us confidence that even 
when subsampled by a factor of 10, MABEL signal strength is 
sufficient to determine ice- sheet slope on ~90-m length scales. 

IV. Conclusion 

As determined in previous studies, MABEL photon-counting 
lidar data have provided precise knowledge of the along-track 
ground surface [4], [5]. Subsequent releases of the data, which 
include beam-angle calibrations, indicate that the MABEL 
ground signal has < 2-m accuracy and ±14-cm precision for 
the April 20, 2012, flight. While instrument accuracy is larger 
than that of other airborne lidars, it is within the MABEL de- 
sign goals (algorithm development, error analysis, and ATLAS 
model validation), where instrument precision is more critical 
to satellite algorithm development. 

Using the release 8 MABEL data, which include beam- angle 
calibrations, we were able to simultaneously analyze multiple 
beams to determine that the lidar accurately characterizes the 
ice-sheet surface, as demonstrated by comparisons of slope on 
^90-m spatial scales with ATM measurements. By analyzing 
data from both variable (Southern Traverse) and flat (Sum- 
mit Area) terrains, we also demonstrate that MABEL slope- 
determination capabilities remain robust through a wide range 
of localized slope assessments. 

Finally, in preparation for the launch of ICESat-2, we as- 
sessed MABEL local across-track slope determination using 
similar beam spacing and a similar transmit-energy- strength 
relationship between a strong and simulated weak-beam pair. 
While the April 2012 MABEL data densities are, at best, 46% 
of the current expectations for ICESat-2, the slope determi- 
nation using a simulated weak beam based on photon counts 
subsampled by a factor of 4 remained robust. However, at a 
subsampling rate of 10, while the mean slope values for both 
study areas remained relatively unchanged from the full-rate 
MABEL ground signal, the changes in the standard deviation of 
the data suggest that the subsampling of the MABEL photons 
is having an adverse effect on the ground- finding algorithm. 

The results based on MABEL give us confidence in the 
ICESat-2 measurement strategy, with respect to beam geometry 
and the transmit-energy- strength relationship. This measure- 
ment strategy allows for the separation of the effects of local 
slope from true ice-sheet elevation change with just two satellite 
passes of a given area. 
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